Tkachenko modes of vortex lattices in rapidly rotating Bose-Einstein condensates 
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We calculate the in-plane modes of the vortex lattice in a rotating Bose condensate from the 
Thomas-Fermi to the mean-field quantum Hall regimes. The Tkachenko mode frequency goes from 
linear in the wavevector, k, for lattice rotational velocities, Q, much smaller than the lowest sound 
wave frequency in a finite system, to quadratic in k in the opposite limit. The system also sup- 
ports an inertial mode of frequency > 2f2. The calculated frequencies are in good agreement with 
recent observations of Tkachenko modes at JILA, and provide evidence for the decrease in the shear 
modulus of the vortex lattice at rapid rotation. 
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The collective modes of the vortex lattice in a rotat- 
ing superfluid have been of considerable interest since the 
1960's. In a classic series of papers Tkachenko [1] showed 
that the lattice supports an elliptically polarized oscilla- 
tory mode, with the semi-major axis of the ellipse orthog- 
onal to the direction of propagation. Such modes were 
observed in superfluid helium in 1982 [2]. Reference [3] 
reformulated the hydrodynamics of rotating superfluids 
to take into account the elasticity of the vortex lattice (in- 
cluding the normal fluid, dissipation, and line bending - 
Kelvin - oscillations of the vortex lines in three dimen- 
sions) and thus describe the Tkachenko modes; effects 
of the oscillations of the vortex lines at finite tempera- 
ture on the long ranged phase correlations of the super- 
fluid were discussed in [4] . The focus of these papers was 
on superfluid helium, where rotational speeds are always 
much smaller than characteristic phonon frequencies, of 
order 1K/H. 

Atomic Bose condensates, on the other hand, allow 
one to study superfluids over a large range of rotational 
speeds [5 9], from the Thomas- Fermi regime where lat- 
tice rotational frequencies, Q, are small compared with 
phonon frequencies - of order the trapping frequency - to 
the quantum Hall limit of rapid rotations [10-13], where 
a condensate in a harmonic trap flattens to a very weakly 
interacting effectively two dimensional system, and rota- 
tional speeds can well exceed phonon frequencies. One 
may identify three distinct physical regimes: first, the 
"stiff" Thomas-Fermi regime, where Q is small compared 
with the lowest compressional frequencies, sfco, where s 
is the sound velocity, fco ~ l/R is the lowest wavenum- 
ber in the finite geometry, and R is the size of the system 
transverse to the rotation axis. In this regime the sys- 
tem responds to rotation effectively as an incompressible 
fluid. At faster rotation, when O 3> sko, but f2 <C ms 2 , 
where m is the atomic mass, the system is in the "soft" 
Thomas-Fermi regime, where compression of the super- 
fluid becomes important in the response of the lattice. 
Finally, when Q 3> ms 2 , the system enters the "mean 



field" quantum Hall regime, in which the condensation 
is only in lowest Landau orbits. Eventually the vortex 
lattice melts [14,15], and the system enters a strongly 
correlated regime [10-12]. 

Observations of Tkachenko modes in Bose-condensed 
87 Rb were reported recently by Coddington et al. [16] at 
rotation speeds up to 0.975 of the transverse trapping 
frequency, w p ; in this regime, fl increases past sko, an d 
although not yet in the quantum Hall regime, the experi- 
ments show effects on the modes frequencies of very rapid 
rotation. The modes in rotating atomic condensates have 
been the subject of several theoretical investigations, in- 
cluding determination of the Tkachenko modes at slow 
rotation taking full account of the finite geometry [17], 
determination of the fundamental modes in the absence 
of effects of the elastic energy of the lattice [18], and 
studies in the quantum Hall regime [15,19-21]. In this 
Letter we derive the modes of the vortex lattice at gen- 
eral rotation speeds. The present analysis is restricted to 
linearized motion in two dimensions transverse to the ro- 
tation axis, and neglects the normal fluid [22]. Our start- 
ing point is the conservation laws and superfluid accel- 
eration equation governing the system, including the full 
elasticity of the lattice. Taking into account the decrease 
of the shear modulus of the vortex lattice with increas- 
ing f2, one can fully understand the measured Tkachenko 
frequencies. 

We consider a rotating gas of bosons described by a re- 
pulsive short range repulsive interaction, U(r) = gS 3 (r), 
with g = 4ira s /m, where a s is the s-wave scattering 
length, in units in which h = 1. The angular momen- 
tum of the superfluid is carried in vortices, which form 
a triangular lattice rotating as a solid body at angular 
velocity Q. We work in the frame corotating with the 
vortex lattice, and denote the deviations of the vortices 
from their home positions by the continuum displacement 
field, e(r, t). In linear order in the vortex displacements, 
the long wavelength superfluid velocity, v(r, t), can be 
written, following [4], in terms of the long wavelength 
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vortex lattice displacement field, e(r,t), and the phase 
<&(r, i) of the order parameter, as 

u + 2!!x( = V$/m. (1) 

This equation follows from the fact that its curl, 

Vx« = -2S1V'£, (2) 

is the conservation law relating the change in vorticity to 
a compression of the vortex lattice, while its longitudinal 
part is trivially the gradient of a scalar [23]. Equation (2) 
constrains the number of degrees of freedom in two di- 
mensions from five (n, v, e) to four. The time derivative 
of Eq. (1) is the superfluid acceleration equation, 



i • 2D = -V(> - K,n ). 



(3) 



where fx is the chemical potential, and, for a harmonic 
confining trap of frequency uo p in two dimensions, 
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Off 



(4) 



In the frame corotating with the vortex lattice, the chem- 
ical potential \i is related to the phase by 



H{r, t) - V cS 



1 d®(r,t) 



ra 



dt 



(5) 



The local elastic energy density of a triangular lattice 
in two dimensions has the form [3] , 



f(r)=2Ci(V-e) 2 + C 2 



dx 



8ty 

dy 



de x <ky 
dy dx 



(6) 



where C\ is the compressional modulus, and C 2 the shear 
modulus of the vortex lattice. In an incompressible fluid, 
C 2 = raO/8 = —C\. The shear modulus C 2 in fact de- 
creases with increasing O, from fin/ 8 in the incompress- 
ible limit O, eventually reaching, in the mean field quan- 
tum Hall limit, the value [15,21], C 2 ~ (81/807r 4 )ms 2 n. 
In this limit C\ — 0. The falloff for small O/tos 2 has the 
form, 



On / n 
Ci ^ — 1 - 7 — ; 



(7) 



a first estimate [21] is 7 <~ 4. 

The dynamics is specified by the superfluid accelera- 
tion equation (3); the continuity equation, 



dn(r, t) 
dt 



+ V-j(r,t) = 0, 



(8) 



where n is the (smoothed) density, and j — nv is the 
particle current; and conservation of momentum: 



dj_ 

dt 



20 x j }+ VP + nW ff = -a, 



(9) 



where P is the pressure; at zero temperature, VP = 
nV/i, while in equilibrium, VP + nVKff = 0. The elas- 
tic stress, a, is given in terms of the total elastic energy, 
P cl = /d 2 r£(r), by 



a{r, t) 
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= -4CiV(V • e) - 2C 2 V 2 e. 



(10) 



Equations (9) and (3), with VP = nV/i, imply 
20 x (e — v) = a j ran. The curl of this equation becomes 



V • (e - v) 



V x a 



C 2 



-V 2 (V x e), 



20nm flnm 
while its divergence, together with (2), yields, 



(11) 



V x e + 20V ■ e = - 



V-cr _ C 2 + 2Ci 
20nm Onm 



V 2 (V-e). (12) 



The structure of the modes follows directly from the lin- 
earized set of equations: (2), (11), (12), and the diver- 
gence of (3). The density oscillations are governed by 



d 2 



n = 2nOV x e, 



(13) 



where the sound speed s is given by ms 2 — dP/dn [24]. 
In the absence of coupling of density oscillations to com- 
pressions of the vortex lattice, this equation is simply 
that of long wavelength phonons in the condensate. Us- 
ing Eq. (11) to eliminate V x e, and ignoring the term of 
order V 4 , we find 



-^ + ^V 2 
dt 2 van 



1 d 2 

n dt z 



(14) 



In the absence of coupling of density oscillations to com- 
pression of the vortex lattice, this equation would de- 
scribe a free Tkachenko mode of the vortex lattice [1,3,4] 
of wavevector k and frequency, lot = (2C 2 /mn) 1 / 2 /c. 

The coupled equations (13) and (14) yield the full spec- 
trum of long wavelength modes. The frequencies, for 
given k, are solutions of the secular equation, 



D{k,u) 



■ LO 



40 2 



s 2 + (C1+C2) k^ 

nra 



+ 



2s 2 C 2lz . 



fc 4 = {oj z -uj){Lo 2 -lo 2 t ) =0. (15) 



For 2s 2 C 2 fc 4 /nm< (40 2 + (s 2 + A(C\ + C 2 )/nm)k 2 ) 2 , as 
is generally the case, the mode frequencies are given by 

4(Ci + C 2 ) 



Luj = 40 2 + 



s 2 + 



nra 



k z 



(16) 



and 
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L0 T 



2C 2 



2 7,4 



S 2 k 



nm (4J1 2 + (s 2 + 4(Ci + C 2 )/nm)k 2 ) ' 



(17) 



The first mode is the standard inertial mode of a ro- 
tating fluid; for £1 <C s 2 k 2 it is a sound wave, while 
for £1 ^> s 2 k 2 , the mode frequencies begin essentially 
at 2f2. The second mode is the observed elliptically po- 
larized Tkachenko mode. See Fig. 1. As Eq. (12) im- 
plies, the longitudinally polarized component is tt/2 out 
of phase with the transversely polarized component; of 
order lot/ '2J1 smaller for the Tkachenko mode and of sim- 
ilar size for the inertial mode [26]. In the stiff limit, the 
Tkachenko frequency, lot — (£l/4m) 1 / 2 fc, is linear in k. 
By contrast, in the very soft limit, the mode frequency, 
lot is quadratic in k at long wavelengths 



LOT 



C 2 
2nm 



1/2 



(18) 



in the quantum Hall regime lot 
[15]. 
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FIG. 1. The inertial and Tkachenko mode frequencies vs. 
wavevector measured in units of 20./ s. The inertial mode 
(upper curve) is in units of 2Q, while the Tkachenko mode 
(lower curve) is in units of (n 3 /ms 2 ) 1/ ' 2 . The modes are cal- 
culated for Ci, C2 in the Thomas- Fermi regime. In the quan- 
tum Hall regime, the Tkachenko modes are softer by a factor 
(9 / V ''){ms 2 / 10 ifi) 
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The very soft Tkachenko mode in the rapidly rotating 
regime leads to infrared singular behavior in the vortex 
transverse displacement-displacement correlations at fi- 
nite temperature, and in the order parameter phase cor- 
relations even at zero temperature [15,21]. In a finite 
system the single particle density matrix, (V'Mr ( r '))> 
falls algebraically as (4N v /ir 2 )~' n for \r — r'\ ~ R, where 
rj ~ n 2 y/lON v /9N , N is the total particle number, and 
N v is the total number of vortices present [21]. However, 
dephasing of the condensate only becomes significant as 
N v — ► N, and not necessarily before the vortex lattice 
melts. 
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FIG. 2. Frequency of the lowest Tkachenko mode. The up- 
per curve is for constant shear modulus, C'2, while the lower 
curve includes a decreasing C2, (Eq. (7)) with 7 = 4. The 
data (triangles), from Ref. [16], are multiplied by a factor 
(iV/2.5xl0 6 )- 2/5 to compare with theory, which is calculated 
for N = 2.5 x 10 6 . The inset shows the Tkachenko frequency 
over the entire range of n; the upper (short dashed) curve 
is the mode frequency to lowest order in the wavevector, the 
middle (dashed) curve includes the full frequency dependence 
at constant C2, and the lower curve includes the decreasing 
C2 as well. 

To compare Eq. (17) with the experiment of Ref. [16], 
we extract the effective value of k in the lowest Tkachenko 
mode from the numerical result of Ref. [17] for the funda- 
mental frequency: lot ^ lA3aCl/R = 2.72(il/m) 1 / 2 /i?, 
where a — (27r) 1 / 2 /3 1/ ' 4 (mri) 1 / 2 is the lattice constant, 
and R the Thomas-Fermi radius in the transverse direc- 
tion. Comparison with the slow rotation Tkachenko fre- 
quency, (il/4m) 1 / 2 A; indicates that the effective wavevec- 
tor of the lowest mode is ko = a/R, where a ~ 5.45, 
the value of k we use in the following. The transverse 
Thomas-Fermi radius of a rotating BEC [25] is 



R 2 = d 2 r(l 



-3/5 



(19) 



where x = f2 2 /cj 2 , d = 1 / ^jmLo p is the transverse oscil- 
lator length, t = l(15Nba s /d)(Lo z /Lo p )} 2 / 5 , and lo z is the 
axial trapping frequency. The sound velocity, evaluated 
in the center of the trap (for simplicity) , is 



ms 



^n(0) = ^r(l - x) 2 / 5 , 



so that 



4/2 



(1-s) 2 / 5 ' 



(20) 



(21) 



Also, (fl/sko) 2 — 2x/(l — x)a 2 : in the present experi- 
ments, Q/sfeo reaches ~ 1.15. Ignoring b ~ 1 here [24], 
as well as the small C\ + C 2 term in Eq. (17), we have, 



,a 2 x l ' 2 {l-xf' b 
'At l-x(l-8/a 2 )' 



(22) 
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The inset in Fig. 2 shows the Tkachenko mode fre- 
quencies as a function of fl/u) p , illustrating the initial 
square root rise and the eventual falloff for ^ u) p . 
The curves are evaluated for the parameters of Ref. [16], 
(u) p ,u) z ) = 27r(8.3,5.2) Hz, at the representative conden- 
sate number, N = 2.5 x 10 6 . Figure 2 proper shows these 
curves in the region measured in [16]. The upper curve is 
calculated from Eq. (22) while the lower includes the de- 
crease (7) in the shear modulus of the vortex lattice, with 
7 = 4. The experiments of [16] are at particle numbers 
- (0.7 - 3) x 10 7 . According to Eq. (22), the frequen- 
cies scale as iV~ 2 / 5 . Thus to facilitate comparison with 
theory, shown for N — 2.5 x 10 6 , we have scaled the in- 
dividual data points down by a factor (N/2.5 x 10 6 ) 2 / 5 , 
which is equivalent to scaling the theory up by the same 
factor. The JILA experiments provide evidence for the 
frequency dependence of the shear modulus of the vortex 
lattice. 
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